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Plates 


Plate  1.  Photograph  of  a  model  of  Hawaii  at  the  City 
of  Refuge  National  Historical  Park,  Hawaii. 
The  color  shading  indicates  the  approximate 
distribution  of  vegetation  and  lava  flows. 

Plate  2.  Photograph  of  the  windswept  trees  southwest 
of  Waimea. 

Plate  3.  Photograph  of  the  Kohala  cloud  from  a  loca- 
tion southwest  of  Waimea. 

Plate  4.  Photograph  of  Hawaii  taken  from  Apollo  9. 
The  NASA  caption  reads:   "Here  the  island 
of  Hawaii  looks  1 i ke  a  dark-visaged  monster 
peering  up  out  of  the  sea  beyond  the  lunar 
module  in  the  foreground.  The  snow-tufted 
peaks  of  Mauna  Loa  and  Mauna  Kea  resemble 
eyes,  and  the  wedge-shaped  cloud  bank  where 
the  trade  wind  was  flowing  over  the  saddle 
between  those  peaks  suggests  a  brow.  There 
are  often  moisture-laden  clouds  on  the  wind- 
ward side  of  a  high  tropical  island  and 
scarcely  any  on  the  leeward  side." 
Arrows  A-A  point  to  the  summit  of  Mauna  Kea; 
arrows  B-B  point  to  the  summit  of  Mauna  Loa. 
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A  THREE-DIMENSIONAL  SIMULATION  OE  WINDS 
AND  NON-PRECIPITATING  OROGRAPHIC  CLOUDS  OVER  HAWAII 


Everett  C.  Nickerson 

and 
Elemer  L.  Magaziner 


A  mesoscale  model  is  presented,  designed  with  the  capability 
of  simulating  moist,  non-precipitating,  three-dimensional  flow 
over  mountainous  terrain.   The  model  makes  use  of  a  transformed 
sigma  coordinate  system  to  assure  adequate  resolution  in  the  bound- 
ary layer.   A  preliminary  evaluation  of  the  performance  of  the 
model  has  been  obtained  from  simulations  of  airflow  and  orograph- 
ically  induced  cloud  cover  over  the  island  of  Hawaii. 


1.   INTRODUCTION 

Numerical  models  on  an  operational  basis  are  now  capable  of  reasonably 
accurate  48-hour  forecasts  of  synoptic  scale  disturbances.  However,  those 
significant  weather  events  that  greatly  affect  our  daily  lives,  such  as  heavy 
precipitation  and  damaging  winds,  are  usually  on  a  much  smaller  scale  and 
therefore  not  directly  resolvable  by  the  large  scale  models.  Detailed  fore- 
casts for  specific  areas  make  use  of  information  contained  in  the  fields  of 
large  scale  convergence,  flow  patterns,  moisture  supply,  and  atmospheric 
stability  as  well  as  local  climatological  information.  That  is,  synoptic 
scale  predictions  of  wind,  temperature,  moisture,  and  pressure  fields  do  not 
in  themselves  indicate  the  occurrence  of  severe  weather  in  a  given  locality, 
but  their  structure  may  indicate  to  an  experienced  person  the  likelihood  of 
such  occurrences. 

If  you  examine  the  status  of  current  modeling  efforts  on  the  scale  of 
individual  convective  clouds,  you  find  that  three-dimensional  time  dependent 
simulations  are  now  possible.  Many  microphysical  questions  need  to  be  re- 
solved and  sub-grid  scale  closure  poses  a  difficult  challenge;  nevertheless, 
the  numerical  simulation  of  isolated  convective  clouds  is  technically  feas- 
ible. 

The  fact  that  significant  progress  has  been  made  in  the  numerical  model- 
ing of  these  two  dissimilar  scales  of  motion  is  not  fortuitous  but  is  rather 
a  manifestation  of  the  energetics  that  characterize  atmospheric  motions. 
Current  evidence  suggests  the  existence  of  a  bimodal  distribution  in  the 
atmospheric  kinetic  energy  spectrum,  with  peaks  roughly  at  4  days  and  at  1 
minute(see  Fiedler  and  Panofsky,  1970).  The  first  peak  is  associated  with 
synoptic  scale  baroclinic  disturbances  and  the  second  with  convective  activity, 


Between  these  two  peaks  lies  a  broad  minimum  or  "mesoscale  gap."  Strong 
controls  on  mesoscale  systems  are  therefore  exerted  by  the  synoptic  scale. 
These  controls,  manifesting  themselves  as  time-dependent  lateral  boundary 
conditions  for  mesoscale  models,  also  provide  the  forcing  that  leads  to 
energy  conversions  on  smaller  scales.  The  specification  of  sub-grid  scale 
processes  is  even  more  difficult  for  mesoscale  than  for  synoptic  scale  mo- 
tions since  the  percentage  of  area  covered  by  active  updrafts  may  not  be  as 
small  as  is  presumed  in  some  parameterization  schemes.  A  mesoscale  model 
must  therefore  provide  realistically  for  larger  scale  forcing  as  well  as  for 
sub-grid  scale  convective  transports  of  a  rather  general  nature. 

Significant  progress  in  the  modeling  of  mesoscale  disturbances  is  essen- 
tial for  more  accurate  forecasts  of  the  damaging  aspects  of  severe  storms. 
Mesoscale  models  are  also  needed  for  the  planning  and  evaluation  of  future 
attempts  to  modify  severe  storms,  since  fiscal  constraints  would  prohibit  the 
gathering  of  statistically  significant  data. 

Cumulus  parameterization  theory  is  being  developed  on  a  very  sophisti- 
cated level  for  general  circulation  models,  and  yet  the  local  recycling  of 
mass  and  momentum  by  individual  clouds  is  still  not  well  understood.  Despite 
the  complexity  of  the  problem,  significant  progress  in  developing  a  mesoscale 
model  can  still  be  made  by  formulating  the  basic  equations  to  treat  stable 
precipitation  and  other  important  physical  processes.  The  details  of  a 
three-dimensional  terrain  model  designed  to  simulate  three-dimensional  air- 
flow over  complex  terrain  are  presented  in  sections  2  to  10  of  this  report. 
Preliminary  results  using  the  island  of  Hawaii  as  a  test  case  are  contained 
in  sections  11  to  14. 


2.  BASIC  EQUATIONS 

In  the  presence  of  favorable  synoptic  scale  conditions,  the  vertical 
fluxes  of  heat,  momentum,  and  moisture  within  the  atmospheric  surface  layer 
can  exert  a  strong  influence  on  the  location,  duration,  and  intensity  of 
severe  convective  storms.  Satellite  pictures  have  clearly  shown  that  pre- 
existing cloud  cover  can  effectively  suppress  surface  heating  over  a  wide 
area  and  thereby  prevent  the  development  of  convective  activity  (Weiss  and 
Purdom,  1974).  Moreover,  variations  in  terrain,  coupled  with  unevenness  in 
surface  heating,  can  produce  slope  winds,  mountain-valley  winds,  and  sea 
breeze  circulation  systems.  A  constant  flux  layer  100  m  or  more  in  depth  is 
not  at  all  uncommon  during  periods  of  strong  insolation;  however,  the  depth 
of  that  layer  may  decrease  by  an  order  of  magnitude  during  nighttime  hours. 

In  order  to  provide  for  increased  vertical  resolution  at  the  lower 
boundary  and  still  to  retain  the  advantages  of  a  uniform  grid,  the  following 
coordinate  transformation  was  incorporated  into  the  model: 


a  =  (4v-v4)/3  (1) 


where  a  is  the  conventional  normalized  pressure  coordinate, 
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The  basic  equations  for  the  Nu  coordinate  system  are  as  follows: 
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4  =  -  Cp6(l+0.61qv) 
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(8) 


If  Ps  is  the  surface  pressure  and  Pj  is  the  pressure  at  the  top  of  the 
model  ,  then 


,  =  ps  -  pt 


(9) 


and 


a  =  (P-PT)/W 


(10) 
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We  also  have  made  the  following  definitions  (these  supplement  those  in 

Appendix  A) : 


U   =  TTU 

V  =  ttv 
o'   =  da/dv 
v  =  dv/dt 
P  =  (P/P0)K 

♦  =  gz 

T*  =  T(l+0.61  q  ) 

S  =  ir(ln  (T/p)  +  Lqv/CpT) 

W     =   Tr(qv  +  qcw) 

L     =  597.3  -  0.566  (T  -  273.16) 


(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 


Here  qv  and  qcw  represent  the  mixing  ratio  for  water  vapor  and  cloud 
water,  respectively.  For  all  calculations  reported  on  in  this  paper,  p0  is 
1013  mb,  and  pT  is  zero. 

Temperature  is  not  explicitly  predicted  by  the  model,  but  must  be  ex- 
tracted from  the  predicted  value  of  the  entropy  variable,  S.  If  the  air  is 
saturated,  the  mixing  ratio  is  a  known  function  of  temperature,  and  (18)  be- 
comes a  transcendental  equation  for  the  temperature,  Ts  corresponding  to  sat- 
uration with  respect  to  liquid  water.  If  on  the  other  hand  the  air  is  unsat- 
urated, the  mixing  ratio  in  (18)  is  replaced  by  W/tt  and  the  temperature  is 
solved  for  directly,  using  the  temperature  at  the  previous  time  step  to  ac- 
count for  the  temperature  dependence  of  the  latent  heat.  In  the  present  for- 
mulation of  the  model  no  conversion  of  cloud  water  to  precipitation  is  in- 
cluded, and  saturation  with  respect  to  liquid  water  is  required  for  any  cloud 
water  to  be  present. 

In  order  to  determine  the  temperature  and  moisture  variables,  we  first 
compute  Ts  and  qvs,  the  saturation  temperature  and  its  corresponding  mixing 
ratio  (the  Newton-Raphson  technique  is  currently  being  used).  Temperature 
and  mixing  ratio  are  then  arrived  at  in  accordance  with  the  following  cri- 
teria: 


qv  "  qvs 

q   =  W/tt 
xw 

T  =  T 


-  q 


vs 


W  >  uq 


VS 


qv  =  W/tt 

q   =  0  )  W  ■;  Trq 

TW  ^VS 

T  =  T 

uns 

The  saturation  vapor  pressure  with  respect  to  water,  es ,  used  to  compute 
qvs  is  taken  from  Murray  (1967), 

e  =  6.11  exp  [17.27(T-273.16)/(T-35.86)1  .  (21) 


3.  LOWER  BOUNDARY  CONDITIONS 

A  constant  flux  layer  is  presumed  to  exist  between  the  lower  boundary 
(with  roughness  length  z0) ,  and  the  first  grid  point  above  the  surface, 
z  =  h.  The  friction  velocity  and  surface  fluxes  of  sensible  and  latent  heat 
are  computed  using  the  Bussinger-Dyer  surface  layer  formulation  (see  Nicker- 
son  and  Smiley,  1975) . 

^=F(r<!0)  ^ 


(e-e0)  u*    /     v 

■a(r-y-  (23) 


-w'6 


0 


where  L  is  the  Monin-Obukhov  length,  and  where  the  functional  forms  for  F 
and  G  depend  on  the  density  stratification. 

In  order  to  link  the  theoretical  formulation  with  the  discretization 
scheme  used  in  the  model  it  is  convenient  to  define  a  bulk  Richardson  number 
for  the  surface  layer, 

hg(e,  -  e0) 
RH  = ,  (24) 

B        VSi 


where  the  subscript  h  denotes  the  value  at  the  height  h, 


If  we  now  set 


L  = 


-u|  e. 


(25) 


kg  w"95 


(25)  may  be  rewritten  in  the  form 


F(z/L)  =  Gz/L  -  k  F2R.  =  0  . 


(26) 


For  a  given  value  of  bulk  Richardson  number  and  surface  roughness,  (26) 
becomes  a  transcendental  equation  for  z/L,  where  z  is  now  h.  After  solving 
(26)  for  h/L,  we  then  have 


u*  =  uh/F 


(27) 


we5  =  uh(e0  -  eh)/FG 


(28) 


w'qo  =  uh(q0  -  qh)/FG 


(29) 


For  unstable  conditions  (i.e.,  e0  -  Q^  >   0),  the  functions  F  and  G  have 
the  form 


kF  =  In 


(c-lKco+D" 

-(ctD(co-i'). 


+  2  tan"1  c  -2  tan"1  5, 


(30) 


kG  =  R  In 


(n2-D(nH) 
L(n2+D(n§-l)J 


(31) 


where 


?     =    (1   -  yz/L)1^      ,  (32) 


Co  ■  (1  -  YZ0/L)l/,f    ,  (33) 


n     =    (1   -  y"z/L)1/t+    ,  (34) 


no  =   (1   -  Y"z0/L)lA    ,  (35) 


For  mildly  stable  conditions  (i.e.,  z/L  <  1),  the  functions  F  and  G  have 
the  form 

kF  =  In   (z/z0)  +  Bz/L  ,  (36) 

kG  =  R  In   (z/z0)  +  6z/L  .  (37) 

For  z/L  >  1,  the  functions  F  and  G  are  obtained  from  profiles  given  by  Webb 
(1970). 

kF  =  3  In   (z/L)  +  In   (z/z0)  +  e  (38) 

kG  =  (1  +  3  -  R)  In   (z/L)  +  R  In   (z/z0)  +  B  .       (39) 

At  the  present  stage  of  model  development  the  surface  temperature  and 
mixing  ratio  are  specified  initially  and  not  allowed  to  change  during  the 
course  of  the  model  run. 


4.   PLANETARY  BOUNDARY  LAYER 

A  boundary  layer  is  presumed  to  exist  between  zg,  the  first  grid  point 
above  the  surface  and  z/\,   some  specified  height  where  the  winds  become  de- 
coupled from  surface  layer  effects.  At  the  present  time  z/\  is  set  equal  to 
one  kilometer.  Exchange  coefficients  for  momentum,  K(u),  and  for  the  temp- 
erature and  moisture,  K(T),  between  z/\  and  Z3  are  computed  in  accordance 
with  the  profile  given  by  O'Brien  (1970). 


K  =  KA  +  [(z-za)/(az)2J  [kb  -  KA  +  (z-zB)(Kg+2  (KR-KA)/Az)]    (40) 

The  prime  denotes  a  derivative  with  respect  to  z,  and  az  =  (z/\  -  zg). 
Both  K(u)  and  K(T)  are  set  equal  to  zero  for  heights  greater  than  z/\.  More- 
over, horizontal  friction  is  not  explicitly  included  in  the  model  except  at 
the  top  and  on  the  lateral  boundaries  (see  section  10). 

The  frictional  terms  in  (2),  (3),  (4),  and  (5)  are  written  as 

Fu  "  A  h    (a«u)  &)  (41) 


K   =  A  It:  (ak(U)  §£  )  (42) 


V      3v 


Fs  ■ A  k  (AK(T)  #)  (43) 

Fw  ■ A  h  (AK(T)  Hf)    •  (44) 


where  A  =  ar^?   • 


At  the  lower  boundary  (v  =  a   =  1), 


AK(u)  |ij-=  tt  U2  cos  a  (45) 


AK(u)  |^=  tt  U2  sin  a  (46) 


AK(T)||=n(V+^e)  (46) 


3W 


AK(T)^=ttQ6  (48) 


where  tan  a  =  V^/U^,  the  ratio  of  the  two  wind  components  at  the  first  grid 
point  above  the  surface,  and  where 


•"Qs  =  (e0-eh)  /uh  +  vh  /FG  <49) 

7rQe  =  (q0-qh)  /U2  +  V2  /FG  (50) 


5.  VERTICAL  GRID 

The  vertical  grid  consists  of  15  equally  spaced  levels  in  the  Nu  coor- 
dinate system.  With  reference  to  Figure  1,  the  variable  v  is  defined  at  the 
circled  levels,  and  all  other  variables  are  defined  at  the  crosses. 

Vertical  distances  in  the  Nu  coordinate  system  range  from  zero  to  unity, 
just  as  in  the  sigma  system.  Furthermore,  equation  (1)  satisfies  the  condi- 
tion noted  by  De  Rivas  (1972),  that  da/dv  must  be  finite  over  the  entire  do- 
main and  must  be  equal  to  zero  at  v  =  a  =  1  in  order  to  assure  second  order 
accuracy  in  the  discretization  scheme.  The  vertial  grid  corresponding  to  the 
15  levels  in  the  model  is  shown  in  Table  1,  together  with  the  corresponding 
sigma  levels.  The  first  grid  point  is  located  approximately  18  meters  above 
the  lower  boundary. 


Table.  7.     CompaJvUon  ofa  the.  Two  Cooh.cU.nate.  SyAtemi'   UeAticat  GiicU 


0.0000  0.0000 

0.0333  0.0444 

0.1000  0.1333 

0.1667  0.2220 

0.2333  0.3101 

0.3000  0.3973 

0.3667  0.4829 

0.4333  0.5660 

0.5000  0.6458 

0.5667  0.7212 

0.6333  0.7908 

0.7000  0.8533 

0.7667  0.9071 

0.8333  0.9504 

0.9000  0.9813 

0.9667  0.9978 

1.0000  1.0000 


6.   HORIZONTAL  GRID 

A  staggered  grid  (fig.  2)  is  used  to  reduce  the  truncation  errors  on  the 
level  Nu  surfaces  (see  Anthes  and  Warner,  1974).  The  variables  u  and  v  are 
defined  at  the  crosses,  while  all  thermodynamic  variables  are  defined  at  the 
dots.  The  variable  v  is  not  defined  on  this  horizontal  level  but  at  dis- 
tances Av/2  above  or  below  this  level  and  at  positions  corresponding  to  the 
dots.  For  the  runs  reported  on  in  this  paper,  Ax  and  Ay  have  both  been  set 
equal  to  10  km.  With  26  grid  points  in  each  direction,  this  corresponds  to 
an  area  250  km  on  a  side. 


wmmmmmmm. 

Grid  Level 

Condition  on  Nu 
i/=0 

X 

1 

o 

X 

2 

o 

X 

3 

o 

X 

4 

• 

• 

• 

X 

13 

o 

X 

14 

o 

Aj/=1/15 

X 

15 

v=\ 

i-1,j+1 


i-1. 


I-1.J-1 


Ax 


i,j  +  1     i  +  1.j  +  1 


i,j        i+1, 


i.j-1    i  +  1.j-l 


,AV 


FiguAe.  I.      VeAllcal  qfuA  used  In  the.         VIquJkl  1.     Horizontal,  qhld  u^ed  In 
model  calculation*.  the.  model  calculations. 
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7.  TIME  DIFFERENCING 

Centered  differences  are  used  to  represent  the  time  derivative,  except 
that  every   fifth  time-step  a  Time  and  Space  Uncentered  Matsuno  procedure  is 
introduced  (see  Arakawa  and  Mintz,  1974).  If,  for  example,  we  are  dealing 
with  the  equation 


oT'F-  <51> 


then  the  TASU  -  Matsuno  scheme  is  as  follows 


A*(n+1)   An      n 
A    At     'FC<A' 

A"*'  -  A"  =  F  (A*(n*D, 

At       urx      ' 

A*(n+2)  .  fln+l  _     . 

At         cv    '  Forward       Centered 


fln+2   fln+l        *fn+:M 

A  A  + =  Fni  (A  ^n+6>)  Backward      Lower-Left 

At 


Time 

Space 

Forward 

Centered 

Backward 

Upper-Right 

8.  SPACE  DIFFERENCING 

The  horizontal  grid  and  the  discretization  scheme  that  follows  are 
based  on  Anthes  and  Warner  (1974).  The  staggered  grid  has  a  smaller  trunca- 
tion error  than  a  non-staggered  grid.  Perhaps  more  important,  however,  is 
the  fact  that  boundary  values  of  the  velocity  components  appear  only  in  the 
flux  terms  of  the  momentum  equations  and  not  in  calculations  using  the  con- 
tinuity equation. 

The  finite  difference  analogs  of  equations  (2)  to  (8)  are 


(w)  =   Di(U)  +  D,(U)  +  D3(U)  +  fV  +  D6(U)  (52) 
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(if)  =  °i(V)  +  D^y)   +  D3(V)   -  fU  +  D7(V) 


(53) 


I)  =  D2(S)   +  DJS)  +  D5(S) 

/      i+^j+^k 


(54) 


(t-£)  =  D2(W)   +  DJW)  +  D5(W) 


(55) 


/     \  15 

(|f)  =  -  Av      Z      a'[D8(U)   +  D9(V) 


(56) 


(vaMi+v.   i+i,  ,^  =  -  Av      £    a'     f£  +  D8(U)   +  D9(V) 
iTisJTiji\ — a  k*=l         *- 


3tt 


(57) 


*H%,j-tfs,k  *    ♦l-Hs,j-HSfk+l  +  CP 


[ek+,(l+0.6lVk+!2(fk+1-Pk)i+^+J(58) 


where  the  following  definitions  have  been  made: 


Di(x)   =  - 


4a 


k     [<U1+1J^1  J>(X1*jV  Xi+k,3-k] 


-   W   L(Vi,J+l+    Vi,j)(ri^,j^a  +   V2,W 

i-^J-^J 


-    (V.     .    +   V.     .    ,)(X.,,      •    ,     +    A 


j_  r 

a'Av      |_ 


D2(X)    =  ^777     |(*a'A)k+Jg  -    ( 


va 


A*] 
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D3(X) 


(Av)2 


(Xk+1  -  Ak)(AK(u))         -   (Ak  -  X^flWu)^. 


D^(A)    =   "  4     [<Ui+l,j+l   +   U1+1J«XW/2.W  +  V^W 
■(U1.W   +  U1.J)(X14%.W+^H|J4I|)] 


(Vi  +  l,j  +  l+Vi,j  +  l^Ai^,j  +  3/2^i+,2,j+,2) 


4Ay 


D*(u)  =  2Z7   ['n^j*  '  ¥1-W*  +  '1-Mk.JJk  "  "i-H,3-H 


-  ( 


^i-^j-*] 


D7(V)   = 


2Ay 


p+Jg.j+Js   "    "l+^j-Jg  +   ^-JS.J^   "    *1-%,J-%"J 


"  W    [(7r*)1^.J+Js  -  (*Wj-%  +  ^l-W*  "  ^i-W-* 


D8(U) 


2ax 


Ji+ij+i-uij+i  +  ui+i,j-ui,j 
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Mv>-^  [vi+i.j+i  "  VM.J.+ vi.j«  "  vu] 


U     =     TTU 


V     =     TTV 


\+k      2  ^Xk  +  Xk+1^ 


Xi+ia,J+%  =  4  (xi+l,j+l  +  Xi+l,j  +  Xi,j+1    +  Xi,j} 


C  =  $  -  RT*cnr/P 


The  vertical  velocity,  dz/dt,  is  also  calculated  in  order  to  assist  in 
the  interpretation  of  the  results.  This  is  especially  important  for  computa- 
tions above  steep  mountain  slopes  where  the  usual  approximation  w  =  -  u/pg 
can  even  give  the  wrong  sign.  The  vertical  velocity  is  given  by 


w 


=  1 
i+%,j+%,k   g 


,  (<|,i+3/2,j+is_"  ^yWk  - 

2ax         i+^j+^k 


,  (*i+%,j+3/2  '  *i+3s,j-3g)k  - 

2Ay  i+^j+^k 


■(ttI^^W,^] 
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9.   SMOOTHING  OPERATOR 

In  order  to  suppress  2ax  noise  in  the  flow  fields,  a  smoother  (see 
Shapiro,  1970)  is  applied  to  the  predicted  fields  of  ttu,  and  ttv. 

A  smooth  value  of  a  variable,  A,  is  defined  by: 


*-F<X1+lJ+l  +Ai-l,j+l+Ai+l,j-l+Xi-l,j-l+4>i,j) 


10.  BOUNDARY  CONDITIONS 

At  the  upper  boundary  (a  =  v  =  0),  v  is  set  to  zero.  This  is  a  reason- 
able condition  when  Pj  =  0,  but  may  lead  to  difficulties  when  Pj  is  greater 
than  100  mb.  Quantities  such  as  wind  and  temperature  are  set  equal  to  values 
at  the  grid  point  immediately  below  the  top  of  the  model.  Actual  values 
chosen  are  not  of  great  importance,  however,  since  the  equations  are  written 
in  flux  form  (i.e.  v-rru)  and  the  terms  containing  boundary  values  vanish  on 
the  upper  boundary  due  to  the  condition  on  v. 

The  thermodynamic  variables  are  specified  at  all  boundaries.  Winds  are 
specified  on  inflow  boundaries  and  extrapolated  from  interior  values  on  out- 
flow boundaries.  In  addition,  a  viscous  term  is  added  to  the  predictive 
equations  for  U  and  V  on  the  lateral  and  top  boundaries.  This  is  applied 
only  at  the  first  two  interior  grid  points  and  is  a  simple  Laplacian  operator 
evaluated  at  the  (n-1)  time  level. 


11.  DISCUSSION 

Model  validation  is  an  important  part  of  any  model  development  program. 
Before  a  model  can  be  accepted  as  a  useful  operational  or  research  tool,  it 
must  first  be  demonstrated  that  the  model  is  capable  of  simulating  signifi- 
cant features  of  the  prototype  data.  The  selection  of  a  validation  site  is 
an  especially  important  step  in  the  evaluation  of  a  complex  meteorological 
model,  because  a  site  with  well-defined  patterns  in  the  winds,  precipitation, 
and  cloud  cover  can  provide  many  opportunities  for  comparison  between  model 
results  and  prototype  data.  The  more  features  that  a  model  is  able  to  repro- 
duce, the  greater  the  confidence  that  can  be  placed  in  calculations  for  loca- 
tions lacking  a  good  data  base.  Moreover,  in  view  of  the  numerous  options 
available  to  the  numerical  modeler,  such  a  site  can  also  assist  in  the  formu- 
lation and  specification  of  computational  boundary  conditions,  smoothing 
operators,  and  sub-gridscale  parameterizations. 

Within  this  context,  the  island  of  Hawaii  (Plate  1,  Figs.  3-4)  appears 
to  be  an  ideal  prototype  for  a  mesoscale  model  development  program.  The 
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relatively  smooth  topography,  extending  from  sea  level  to  4.2  km,  is  charac- 
terized by  barren  lava  fields,  dry  grasslands,  rolling  fields  of  sugar  cane, 
and  tropical  rain  forests.  Considerable  data  have  been  gathered  on  rainfall 
and  other  meteorological  parameters  during  Project  Shower  (Mordy  et  al . , 
1957),  and  also  during  the  Warm  Rain  Project  (Lavoie,  1966,  and  Lavoie  et 
al.,  1967).  Part  of  a  concise  introduction  to  the  meteorology  of  the  island 
(NOAA,  1972)  is  reproduced  in  Appendix  C. 

The  importance  of  the  twin  peaks  Mauna  Loa  and  Mauna  Kea  to  the  island 
meteorology  cannot  be  overemphasized.  Leeward  convergence  of  the  anti- 
cyclonic  flow  south  of  Mauna  Loa  and  the  cyclonic  flow  north  of  Mauna  Kea 
results  in  a  suppression  of  the  vortex  streets  commonly  observed  in  the  cloud 
patterns  downwind  of  tall  islands.  This  zone  of  horizontal  convergence  (see 
Fig.  5,  taken  from  Patzert,  1970)  apparently  interacts  with  the  sea  breeze  to 
produce  a  secondary  rainfall  maximum  on  the  leeward  (Kona)  coast  (Chopra, 
1973).  The  Kona  rainfall  maximum  and  the  sharp  demarcation  between  wet  and 
dry  areas  are  vividly  illustrated  in  Fig.  6,  which  shows  Lavoie' s  (1966) 


VigvJie  4.     TenJuzin  peupec- 
tive  oft  Hawati.     Locations 
ol  individual.  gnJA  points 
coincide  with  the  inteh.- 
i>ectioni>  ofa  north-south  and 
eaAt-weAt  lines.     The  view 
ij>  ^rom  northeast   (50°) 
looking  southwest. 
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summer  rainfall  distribution.  Despite  the  occurrence  of  different  wind 
regimes  associated  with  winter  storm  systems,  the  locations  of  the  precipi 
tation  maxima  and  minima  in  Fig.  6  also  coincide  with  those  of  the  annual 
pattern. 
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12.   INITIAL  CONDITIONS 

During  the  fully  developed  and  undisturbed  trade  wind  regime,  the  air  in 
the  lowest  few  kilometers  of  the  Eastern  Pacific  high  has  ample  time  to 
adjust  to  conditions  at  the  ocean  surface  before  encountering  land.  Under 
those  conditions,  topographical  features  or  other  sources  of  horizontal 
variability  outside  the  computational  domain  ought  to  have  little  influence 
on  the  island  weather.  Representative  inflow  boundary  data  could  therefore 
be  obtained  from  an  undisturbed  trade  wind  sounding.  The  sounding  station  at 
Lihue,  Kauai  would  probably  provide  the  most  representative  data  on  a  routine 
basis. 

Initial  conditions  of  wind,  temperature,  and  mixing  ratio  for  the  numer- 
ical calculations  are  shown  in  Figs.  7,  8,  and  9.  The  initial  wind  field 
contains  easterlies  at  low  levels  and  westerlies  aloft.  The  temperature  at 
the  lower  boundary  remains  constant  in  time,  and  decreases  linearly  with 
altitude  from  its  sea  level  value  of  299°K  (see  Fig.  10). 

The  temperature  at  the  lower  boundary  has  not  been  related  to  land  use 
or  vegetation;  however,  limited  use  of  that  information  has  been  made  in 
specifying  the  roughness  length  and  mixing  ratio  at  the  surface.  Vegetation 
and  soil  conditions  were  estimated  from  U.S.G.S.  woodland  print  topographic 
charts.  Roughness  lengths  were  then  obtained  from  Priestley  (1959),  Sellers 
(1966),  and  Sutton  (1953).  Terrain  heights  for  individual  grid  points  to- 
gether with  the  corresponding  surface  parameters  are  listed  in  Appendix  B. 


13.  MODEL  RESULTS 

All  model  results  reported  in  this  section  were  calculated  in  a  terrain- 
following  coordinate  system.  However,  it  is  often  more  convenient  to  display 
the  data  using  height  as  the  vertical  coordinate.  Since  the  height  of  each 
grid  point  is  known,  the  value  of  a  guantity  at  some  intermediate  location  is 
calculated  using  linear  interpolation  between  two  computational  levels. 

Horizontal  windfields  corresponding  to  elevations  of  100,  1000,  1500, 
2500,  and  4000  meters  above  sea  level  are  shown  in  Figs.  11-15.  The  low 
level  convergence  zone  on  the  left  side  of  Fig.  11  is  of  particular  interest, 
and  an  expanded  view  of  the  convergence  zone  is  shown  in  Fig.  16.  Although 
the  computed  winds  do  not  contain  all  the  details  of  the  composite  wind  field 
shown  in  Fig.  5,  the  model  did  succeed  in  reproducing  some  of  the  primary 
features  including  the  leeward  convergence  zone. 

The  numerical  computations  indicate  that  the  effects  of  the  island  on 
the  airflow  may  be  felt  50  to  100  km  from  land.  A  more  accurate  determina- 
tion of  the  island  influence  region  would  reguire  the  calculations  to  be 
performed  over  a  larger  horizontal  domain. 
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For  the  particular  set  of  initial  conditions  used,  Fig.  17  indicates 
that  the  strongest  surface  winds  are  located  near  the  summit  of  Mauna  Kea. 
At  lower  elevations,  the  area  of  strongest  surface  winds  is  along  the  north 
shore  area  and  in  the  saddle  region  between  Mauna  Kea  and  the  Kohala  moun- 
tains (hereafter  called  the  Kohala  saddle).  The  asymmetry  in  the  wind  field 
is  very  apparent  to  the  casual  observer  of  Hawaiian  weather.  With  the  pos- 
sible exception  of  some  promontories,  nowhere  along  the  entire  main  highway 
that  circles  the  island  is  one  more  aware  of  the  presence  of  persistent 
strong  winds  than  in  the  region  just  south  of  Waimea  (the  leeward  edge  of  the 
Kohala  saddle).  Strong  winds  may  occur  at  other  locations  but  the  trees  of 
that  area  vividly  portray  their  struggle  for  survival  in  the  face  of  the 
persistent  and  forceful  northeasterlies  (see  Plate  2). 

Primary  regions  of  convergence  and  divergence  are  quite  apparent  in  the 
horizontal  wind  fields  shown  in  Figs.  11-15.  And  yet  there  must  be  other, 
more  subtle  horizontal  convergence  patterns  responsible  for  the  extreme 
variability  in  the  precipitation  distribution  of  Fig.  6.  Those  terrain- 
induced  regions  of  horizontal  convergence  become  much  more  obvious  when  we 
examine  the  vertical  motion  field  at  an  elevation  of  1000  m  above  sea  level 
as  shown  in  Figs.  18a, b,c.    An  examination  of  Fig.  6  shows  that  all  the 
major  precipitation  maxima  are  well  correlated  with  the  areas  of  rising 
motion  in  Figure  18a.  Areas  of  strong  descent  in  Fig.  18b  also  have  their 
counterparts  in  the  precipitation  minima  of  Fig.  6.  The  Kau  and  South  Kohala 
desert  areas  in  particular  are  under  the  influence  of  sinking  air,  although 
the  model  calculations  seem  to  have  the  North  Kona  rising  motion  intruding 
upon  the  South  Kohala  desert.  The  orographically  forced  upward  motion  on  the 
windward  slopes  is  relatively  shallow  (see  Fig.  19),  but  in  Fig.  20  we  see 
that  the  leeward  convergence  zone  is  associated  with  a  broad  area  of  gradual 
ascent. 

The  correlation  between  the  vertical  motion  field  and  observed  precipi- 
tation does  not  carry  over  to  the  calculated  cloud  cover  distribution  nearly 
as  well.  A  vertical  cross  section  through  the  Mauna  Loa-Mauna  Kea  saddle 
cloud  (Fig.  21)  shows  a  cloud  only  half  the  thickness  of  that  observed  by 
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ViguAe  IS.     a.   Voiijtive  ventrical 
velocity  {ield  1000  m  above.  i>ea 
level  a^tex  5015  6  ofi  model  time.. 
Shaded  a/ieaA  indicate  valaei  in 
exceAA  0^  5  am/6,     b.   Negative 
vertical  velocity  iield  1000  m 
above  tea  level  OL^tQA  5015  4  o{> 
model,  time..     Shaded  aAeat,  indi- 
cjat.2.  6i.nking  motion  in  ex.cei>i>  o{, 
5  cm/6,     c.   Vertical  velocity  at 
1000  m  above  hea  level  afcteh. 
5015  -6  o{  model  twe. 


Mordy  et  al.    (1957).     It  would  probably  be  necessary  to  include  warm  rain 
processes  and  a  convective  parameterization  scheme  in  order  to  simulate  the 
observed  cloud  cover  more  realistically.     The  lack  of  any  cloud  cover  in  the 
Kona  convergence  zone  should  probably  not  cause  any  great  surprise  in  view  of 
the  model's  neglect  of  cumulus  processes.     The  calculated  North  Kohala  cloud 
(Figs.   22-23),  on'  the  other  hand,  appears  deeper  than  the  saddle  cloud  and 
may  be  representative  of  the  actual   Kohala  cloud   (Plate  3). 
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o{  vertical  veloctty  thAough  the 
saddle,  cloud  a^ten.  5015  6  ofi  model 
time..     Shaded  an.ea6  indicate  valuer 
In  exce66  oh  20  cm/6. 


Tlguxe  10.     We6t-ea6t  cn.o66  section 
o&  vertical  velocity  through  the 
Kona  convcn.go.nc2.  zone  a^ten.  5015  6 
0^  model  time.     Shaded  a/iea6  Indi- 
cate valuer  -in  exce66  o£  20  cm/6. 
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Vigu/ie  21.     Cloud  wateA  mixing  latio 
thAough  the  6addle  cloud  cokac6- 
pondlng  to  the  6ame  cao66  6ection 
o6  Vig.    7  9.     Shaded  aAeo6  Indicate 
valuc6  In  exce66  ojJ   7  gm/kgm. 


FiguAe  22.     HoAth-6outh  CA066  6ectlon 
o^  cloud  uxuteA  6houiing  a  pontic n 
o&  the  Kolvxla  cloud  a^tcA  5015  6 
o\  model  tone.     Shaded  oAea6   indi- 
cate value*  in  exce66  OjJ   7  gm/kgm. 
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VIquaz  23.     An  expanded  vim  oft  the 
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Fig.   22.     The  6 haded  axea  Indicate* 
value*  In  exceA*  o{  1  gm/kgm. 
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14.  CONCLUSION 

Terrain-following  coordinate  systems,  such  as  the  one  used  in  the  pres- 
ent calculations,  eliminate  the  need  for  special  boundary  treatments  where 
quasi -horizontal  computational  surfaces  intersect  the  surface  topography. 
However,  one  seldom  gets  something  for  nothing.  In  the  case  of  a  non-uni- 
formly  spaced  terrain-following  system,  the  modeler  must  contend  with  vari- 
able vertical  resolution.  The  vertical  resolution  shown  in  Table  1  is  prob- 
ably sufficient  for  most  applications  in  the  lower  troposphere.  On  the  other 
hand,  upper  tropospheric  or  stratospheric  applications  would  require  a  differ- 
ent coordinate  transformation  in  order  to  match  the  region  of  primary  inter- 
est with  the  region  of  improved  vertical  resolution,  which  occurs  near  the 
lower  boundary  in  the  Nu  system. 

In  spite  of  the  above  caveat,  preliminary  results  indicate  that  the  Nu 
coordinate  system  model  is  capable  of  simulating  not  only  the  significant 
features  of  the  wind  field  (i.e.  the  Kona  convergence  zone),  but  also  the 
smaller  scale  features  that  correspond  well  with  observed  precipitation  pat- 
terns. Apparently  the  smoothing  of  the  pressure  weighted  wind  field  still 
allows  the  model  to  resolve  those  important  smaller  scale  features. 

The  importance  of  the  island  of  Hawaii  as  a  prototype  for  mesoscale 
model  development  cannot  be  overemphasized.  Considerable  data  already  exist, 
and  reasonably  comprehensive  sets  of  initial  and  boundary  data  could  be  ob- 
tained at  relatively  small  additional  cost.  However,  new  interactive  pro- 
grams combining  numerical  modeling  studies  and  field  investigations  would 
derive  considerable  benefit  from  the  remote  sensing  capabilities  of  aircraft 
and  satellites  (see  Plates  1  and  4).  Comprehensive  studies  of  the  spatially 
varying  but  semi -permanent  wind,  cloud,  and  precipitation  regimes  of  Hawaii 
would  (1)  assist  in  the  evaluation  of  mesoscale  models,  and  (2)  provide  the 
physical  insight  required  for  improved  parameterization  of  sub-grid  scale 
transport  processes. 
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APPENDIX  A. 
BASIC  NOMENCLATURE 


w 

G 

K(u) 

K(T) 

L 
L 
P 
P 
P 


R 

R 

R. 
*l 
S 

T 

T. 


uns 


Unspecified  variable  at  time  level  n 

Specific  heat  at  constant  pressure 

Normalized  surface  layer  wind  speed 

Friction  term  in  the  entropy  equation 

Friction  term  in  the  x  equation  of  motion 

Friction  term  in  the  y  equation  of  motion 

Fraction  term  in  the  moisture  equation 

Normalized  surface  layer  heat  flux 

Exchange  coefficient  for  momentum 

Exchange  coefficient  for  thermodynamic 
variables 


Monin-Obukhov  length 
Latent  heat  of  vaporization 
Pressure 

Normalized  pressure 
Reference  pressure 


uie, 


-m 


kg  w'8g 


Pressure  at  the  ground  (v  =  1) 

Pressure  at  the  top  of  the  atmosphere 

(v  =  0) 

Moisture  flux 

Sensible  heat  flux 

Gas  constant   [=  287.04] 

Boundary  layer  constant   [=    .74) 

Bulk  Richardson  number 

Entropy  variable     [=  it(ln   (T/p)   +  Lq  /CpT; 

Temperature 

Temperature  corresponding  to  saturated 
conditions 

Temperature  corresponding  to  unsaturated 
conditions 

Virtual  temperature  [=  T ( 1+0. 61  q  )] 

Pressure  weighted  wind  component  [=nu] 

Pressure  weighted  wind  component  at  z=h 

Pressure  weighted  wind  component  |=nv] 

Pressure  weighted  wind  component  at  z=h 

Pressure  weighted  water  mixinq  ratio 

[=  n(q  +q   )  ] 
1    VMv  Mcw 


f 

g 

h 
i.j.k 

k 

qcw 

qh 

% 

\ 

\s 

t 

u 

uh 

V 

w 
x,y 

z 

ZA 
ZB 

zo 

Wq6 

w'ej 

a 

B 

Y 

y" 

c 

CO 

n 

no 


Saturation  vapor  pressure  with  respect 
to  water 

Coriolis  parameter  [=  5x10"^ ] 

Acceleration  of  gravity   [=  9.8062] 

Heiqht  of  the  first  grid  point  above  the 
lower  boundary 

Grid  locations  in  the  (x,y,v)  coordinate 
system 

Von  Karman  constant  [=  .35] 

Cloud  water  mixinq  ratio 

Water  vapor  mixing  ratio  at  z=h 

Water  vapor  mixing  ratio  at  z=zQ 

Water  vapor  mixing  ratio 

Saturated  mixing  ratio 

Time 

Velocity  component  in  the  x  direction 

Friction  velocity 

Wind  speed  at  z  =  h 

Velocity  component  in  the  y  direction 

Velocity  component  in  the  z  direction 

Horizontal  coordinates  in  the  Nu  coord- 
inate system 

Vertical  coordinate 

Height  of  the  planetary  boundary  layer 

Heiqht  of  the  surface  layer 

Surface  roughness  length 

Surface  moisture  flux 

Surface  heat  flux 

Surface  wind  direction   [=   tan-1 (Vh/Uh)  ] 

Surface  layer  constant   [=  4.7] 

Surface  layer  constant   (=   15] 

Surface  layer  constant   [=  9] 

Surface  layer  function   [=   (1-yz/L)' 

Surface  layer  function   [=   (l->z:/L)' 

Surface  layer  function   [■   ( 1  - . " z / L ) '    . 

Surface  layer  function   [■   (l-.'z-'L' 
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e  Potential  temperature  [=  T/P] 

6.  Potential  temperature  at  z  =  h 

e0  Potential  temperature  at  z  =  z 

k  Thermodynamic  constant  [=  2/7] 

X  Dummy  variable 

v  Transformed  vertical  coordinate 

v  Time  derivative  of  v;v  =  dv/dt 

tt  Pressure  variable;  tt  =  P_  -  P_ 

o  Vertical  coordinate;  a  =  a   (v)  =  (p-p-r)/iT 

F  Function 

a'  Vertical  derivative  of  o;o'   =  da/dv 

d  Time  derivative  of  o;6  -   da/dt 

<t>  Geopotential ;  4>  =  gz 

10  Vertical  motion  in  the  P  system;  uj=dp/dt 

At  Time  step  [=  10  sec,  except  that  the 
last  time  step  was  15  sec] 

ax,  Ay    Horizontal  grid  length  [=  10  km] 

Az  Boundary  layer  thickness  Az  =  zA  -  zg 

Av  Vertical  grid  length  [=  1/15] 
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APPENDIX  B.  TERRAIN  PARAMETERS 

A  soil  parameter  and  a  vegetation  parameter  have  been  assigned  to  each 
grid  point  at  the  lower  boundary  of  the  model.  The  two  surface  parameters 
are  further  subdivided  and  a  number  is  assigned  to  each  category.  Those  num- 
bers are  then  used  to  estimate  the  surface  roughness  length  and  mixing  ratio 
at  the  individual  grid  points. 


Table.   87.  SoU.  PoAameteJU 

Soil  Type  Coded  Form 

Sea  1 

Dry  2 

Semi  moist  3 

Wet  A 

Lava  5 

Sand  6 


Table  82.      Vegetation  VahjmeteAA 

Vegetation  Type  Coded  Form 

None  1 

Short  Grass  2 

Tal 1  Grass  3 

Shrub  k 

Forest  5 


The  relative  humidity  at  the  surface  was  set  equal  to  1002  except  for 
soil  types  2  and  3,  which  were  assigned  the  values  80"'  and  90'  respectively, 
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Table.  B3.     RoughnzAA  Pcviam&teAA 


Soil    Type  Vegetation   Type  Roughness    Length 

(meters) 

Sea  None  0.0001 

Lava  None  0.01 

Dry  Grass  0.01 

Dry  Shrub  0.2 

Semimoist  Grass  0.05 

Semimoist  Shrub  0.5 

Wet  Grass  0.1 

Wet  Shrub  1.0 

Wet  Forest  3-0 


All  points  on  the  26  by  26  not  listed  below  are  over  water.  Their  sur- 
face parameters  are  coded  (0,1,1),  where  the  first  digit  represent  height 
above  sea  level  in  meters,  the  second  the  soil  parameter,  and  the  third  the 
vegetation  parameter. 
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Table.  B4 .     TeAAain  H(Ugkt6  and  SuA^ace.  ?aAamoJ:&a>    fan.  Hawaii 
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APPENDIX  C.  EXCERPT  FROM  U.S.W.B.  CLIMATOLOGICAL  SUMMARY  (NOAA,  1972) 


The  city  of  Hi lo  is  located  near  the  midpoint  of  the  eastern  shore  of 
the  Island  of  Hawaii.  This  island  is  by  far  the  largest  of  the  Hawaiian 
group,  with  an  area  of  4,038  square  miles,  more  than  twice  that  of  all  the 
other  islands  combined.  Its  topography  is  dominated  by  the  great  volcanic 
masses  of  Mauna  Loa  and  Mauna  Kea,  both  of  which  exceed  13,600  feet  in  ele- 
vation, and  of  Hualalai  (8,271  feet),  the  Kohala  Mountains  (5,480  feet),  and 
Kilauea  (4,090  feet).  In  fact  the  island  consists  entirely  of  the  slopes  of 
these  mountains  and  of  the  broad  saddles  between  them. 

Hawaii  is  diamond-shaped,  about  93  miles  long,  from  north  to  south,  and 
76  miles  wide.  Its  highest  point  is  the  summit  of  Mauna  Kea  at  13,796  feet. 
Mauna  Loa  and  Kilauea,  which  occupy  the  southern  half  of  the  island,  are 
still  active  volcanoes  and  hence  smooth-sloped,  in  contrast  to  the  deeply 
eroded  valleys  that  indent  portions  of  Mauna  Kea  and  the  Kohala  Mountains. 

Hawaii  lies  well  within  the  belt  of  northeasterly  trade  winds  generated 
by  the  semi -permanent  Pacific  high  pressure  cell  to  the  north  and  east.  The 
climate  of  the  island  is  greatly  influenced  by  terrain.  Its  outstanding  fea- 
tures are  the  marked  variations  in  rainfall  with  elevation  and  from  place  to 
place,  the  persistent  northeasterly  trade  winds  in  areas  exposed  to  them,  and 
the  equable  temperatures  from  day  to  day  and  season  to  season  in  localities 
near  sea  level . 

Over  the  island's  windward  slopes,  rainfall  occurs  principally  in  the 
form  of  showers  within  the  ascending  moist  trade  winds.  Mean  annual  rainfall 
increases  from  100  inches  or  more  along  the  coasts  to  a  maximum  of  over  300 
inches  at  elevations  of  2,000  to  3,000  feet,  and  then  declines  to  about  15 
inches  at  the  summits  of  Mauna  Kea  and  Mauna  Loa.  In  general,  leeward 
(southern  and  western)  areas  are  topographically  sheltered  from  the  trades  — 
hence  from  trade-wind  showers  — and  are  therefore  drier;  although  sea 
breezes  created  by  daytime  heating  of  the  land  move  onshore  and  upslope, 
causing  afternoon  and  evening  cloudiness  and  showers.  Where  mountain  slopes 
are  steeper,  mean  annual  rainfall  may  range  from  30  inches  along  the  coast  to 
120  inches  at  elevations  of  2,500  to  3,000  feet.  The  driest  locality  on  the 
island  — and  in  the  State  — with  an  average  annual  rainfall  of  less  than  10 
inches,  is  the  coastal  strip  just  leeward  of  the  southern  portion  of  the 
Kohala  Mountains  and  of  the  saddle  between  the  Kohalas  and  Mauna  Kea. 

These  marked  contrasts  in  rainfall  are  reflected  in  soil  and  vegetation, 
with  frequent  abrupt  transitions  from  lush  tropical  growth  to  near-desert 
conditions,  such  as  occurs  between  Kilauea's  wet  windward  slopes  and  the  Kau 
Desert  just  to  the  south. 

Within  the  city  of  Hi lo  itself,  average  rainfall  varies  from  about  130 
inches  a  year   near  the  shore  to  as  much  as  200  inches  in  mountain  sections. 
The  wettest  part  of  the  island,  with  a  mean  annual  rainfall  exceeding  300 
inches,  lies  about  6  miles  upslope  from  the  city  limits.  Rain  falls  on  about 
280  days  a  year  in  the  Hilo  area. 
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Hawaii's  equable  temperatures  are  associated  with  its  mid-ocean  location 
and  the  small  seasonal  variation  in  the  amount  of  energy  received  from  the 
sun.  At  Hilo,  the  range  in  average  temperature  from  February  and  March,  the 
coldest  months,  to  August,  the  warmest,  is  only  5.2°  F  and  the  average  daily 
range,  15.1°  F.  The  highest  temperature  of  record  at  Hilo  Airport  is  94°  F; 
the  lowest  53°  F.  Greater  variations  occur  in  localities  with  less  rain  and 
cloud,  but  temperatures  in  the  mid-90's  and  low  50's  are  uncommon  anywhere  on 
the  island  near  sea  level. 

The  trade  winds  prevail  throughout  the  year  (although  they  may  be  absent 
for  days  or  even  weeks  at  a  time)  and  profoundly  influence  the  climate.  How- 
ever, the  island's  entire  western  coast  is  sheltered  from  the  trades  by  high 
mountains,  except  that  unusually  strong  trade  winds  may  sweep  through  the 
relatively  low  (2,600-foot)  saddle  between  the  Kohala  Mountains  and  Mauna  Kea 
and  reach  the  areas  to  the  lee.  But  even  places  exposed  to  the  trades  may  be 
affected  by  local  mountain  circulations.  For  example,  the  prevailing  wind  at 
Hilo  Airport  is  not  the  northeasterly  trade,  but  the  southwesterly  wind  that 
drifts  downslope  off  Mauna  Loa  during  the  night  and  early  morning  hours. 

Except  for  heavy  rain,  really  bad  weather  seldom  occurs.  Thunderstorms 
average  only  8  per  year,  and  are  rarely  severe.  During  the  winter,  cold 
fronts  or  the  cyclonic  storms  of  subtropical  origin  (the  so-called  Kona 
storms)  may  bring  blizzards  to  the  upper  slopes  of  Mauna  Loa  and  Mauna  Kea, 
with  snow  extending  at  times  to  9,000  feet  or  below  and  icing  nearer  the 
summit. 

Storms  crossing  the  Pacific  a  thousand  miles  to  the  north,  or  Kona 
storms  closer  by,  may  generate  seas  that  cause  heavy  swell  and  surf  along  the 
northern,  eastern,  and  southwestern  shores  of  the  island. 
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LABOR  AT  DRIES 

The  mission  of  the  Environmental  Research  Laboratories  (ERL)  is  to  conduct  an  integrated  program  of  fundamental 
research,  related  technology  development,  and  services  to  improve  understanding  and  prediction  of  the  geophysical 
environment  comprising  the  oceans  and  inland  waters,  the  lower  and  upper  atmosphere,  the  space  environment  and  the 
Earth.  The  following  participate  in  the  ERL  missions: 


MESA         Marine  EcoSystems  Analysis  Program.  Plans, 
directs,  and  coordinates  the  regional  projects 
of  NOAA  and  other  federal  agencies  to 
assess  the  effect  of  ocean  dumping,  municipal 
and  industrial  waste  discharge,  deep  ocean 
mining,  and  similar  activities  on  marine 
ecosystems. 

OCSEA       Outer  Continental  Shelf  Environmental 
Assessment  Program  Plans,  directs,  and 
coordinates  research  of  federal,  state,  and 
private  institutions  to  assess  the  primary 
environmental  impact  of  developing  petroleum 
and  other  energy  resources  along  the  outer 
continental  shelf  of  the  United  States. 

WM  Weather  Modification  Program  Office.  Plans, 

directs,  and  coordinates  research  within  ERL 
relating  to  precipitation  enhancement  and 
mitigation  of  severe  storms.  Its  National 
Hurricane  and  Experimental  Meteorology 
Laboratory  (NHEML)  studies  hurricane  and 
tropical  cumulus  systems  to  experiment  with 
methods  for  their  beneficial  modification  and 
to  develop  techniques  for  better  forecasting 
of  tropical  weather.  The  Research  Facilities 
Center  (RFC)  maintains  and  operates 
aircraft  and  aircraft  instrumentation  for 
research  programs  of  ERL  and  other  govern- 
ment agencies. 

aoml         Atlantic  Oceanographic  and  Meteorological 
Laboratories.  Studies  the  physical,  chemical, 
and  geological  characteristics  and  processes 
of  the  ocean  waters,  the  sea  floor,  and  the 
atmosphere  above  the  ocean 

PMEL  Pacific  Marine  Environmental  Laboratory. 

Monitors  and  predicts  the  physical  and 
biological  effects  of  man's  activities  on 
Pacific  Coast  estuarine,  coastal,  deep-ocean, 
and  near-shore  marine  environments. 

GLERL        Great  Lakes  Environmental  Research  Labora- 
tory. Studies  hydrology,  waves,  currents,  lake 
levels,  biological  and  chemical  processes, 
and  lake-air  interaction  in  the  Great  Lakes  and 
their  watersheds;  forecasts  lake  ice  conditions 


GFDL      Geophysical  Fluid  Dynamics  Laboratory. 
Studies  the  dynamics  of  geophysical  fluid 
systems  (the  atmosphere,  the  hydrosphere, 
and  the  cryosphere)  through  theoretical 
analysis  and  numerical  simulation  using  power- 
ful, high-speed  digital  computers 

APCL      Atmospheric  Physics  and  Chemistry  Labora- 
tory. Studies  cloud  and  precipitation  physics, 
chemical  and  particulate  composition  of  the 
atmosphere,  atmospheric  electricity,  and 
atmospheric  heat  transfer  with  focus  on 
developing  methods  of  beneficial  weather 
modification 

NSSL       National  Severe  Storms  Laboratory  Studies 
severe-storm  circulation  and  dynamics,  and 
develops  techniques  to  detect  and  predict 
tornadoes,  thunderstorms,  and  squall  lines 

WPL         Wave  Propagation  Laboratory.  Studies  the 
propagation  of  sound  waves  and  electro- 
magnetic waves  at  millimeter,  infrared,  and 
optical  frequencies  to  develop  new  methods 
for  remote  measuring  of  the  geophysical 
environment 

ARL         Air  Resources  Laboratories.  Studies  the 

diffusion,  transport,  and  dissipation  of  atmos- 
pheric pollutants,  develops  methods  of 
predicting  and  controlling  atmospheric  pollu- 
tion; monitors  the  global  physical  environment 
to  detect  climatic  change 

AL  Aeronomy  Laboratory.  Studies  the  physical 

and  chemical  processes  of  the  stratosphere, 
ionosphere,  and  exosphere  of  the  Earth  and 
other  planets,  and  their  effect  on  high-altitude 
meteorological  phenomena 

SEL  Space  Environment  Laboratory  Studies 

solar-terrestnal  physics  (interplanetary,  mag- 
netosphenc.  and  ionospheric),  develops  tech- 
niques for  forecasting  solar  disturbances, 
provides  real-time  monitoring  and  forecasting 
of  the  space  environment 
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